Edge and waveguide THz surface plasmon modes in graphene micro-ribbons 
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Surface plasmon modes supported by graphene ribbon waveguides are studied and classified. The 
properties of both modes with the field concentration within the ribbon area (waveguiding modes) 
and on the edges (edge modes) are discussed. The waveguide and edge modes are shown to be 
separated from each other by a gap in wave numbers. The even-parity hybridized edge mode results 
to be the fundamental electromagnetic mode of the ribbon, possessing also the lowest losses. All 
the plasmonic modes in the ribbons have an optimum frequency, at which the absorption losses are 
minimum, due to competition between the plasmon confinement and the frequency dependence of 
absorption in graphene. 
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The science of graphene is advancing rapidly, discov- 
ering new fundamental physical effects^ and showing 
promising perspectives in several state-of-the-art tech- 
nological applications^. For example, graphene can 
carry both electromagnetic (EM) signals and electric 
currents through the same extremely thin circuitry, 
which could allow electrical interconnects to achieve 
data transmission faster rates. In this regard, it has 
been recently shown that graphene ribbons can operate 
as a broadband radio-frequency mixer at frequencies 
up to 10 GHz inside an integrated circuit^, showing a 
great potential for waveguiding. The two main charac- 
teristics of graphene for this kind of applications are: 
(i) that despite being almost transparent, graphene can 
support surface plasmons (GSPs) in the THz regime^ 
and (ii) the carrier concentration (and thus the conduc- 
tivity) in graphene can be controlled through electro- 
static gating. This last property can, in turn, be used 
to tune the properties of GSPs, opening up a wealth 
of interesting applications in photonics^ s 14^ 

For most functionalities the GSPs must ideally also 
be confined laterally in the graphen e sheet , which is 
known to occur in graphene ribbons^^^^E^EU^ However, 
a characterization of the confined plasmons in this ge- 
ometry is lacking. In this paper we address this prob- 
lem, and study the characteristics of GSP modes in 
graphene ribbons of micrometric widths, in the THz 
regime. We concentrate in the spectral regimes where 
the ribbon supports either a single EM mode or a few 
of them and, in this last case, on the properties of 
strongly localized edge GSP modes (EGSP) that ap- 
pear in the system. The dependence of both waveg- 
uide GSP (WGSP) and EGSP upon relaxation time of 
electrons, permittivity of a substrate and width of the 
ribbon, is also analyzed. 

Let us consider a graphene ribbon of width w (at 




FIG. 1: (Color online) The geometry of the studied system: 
a graphene ribbon of the width w, and conductivity a, 
placed at the interface between two dielectric half-spaces 
with dielectric constants ei and 62- The mode propagates 
along the Oy axis. 



\x\ < 1^/2), placed at the boundary z = between 
two dielectric half-spaces, see Fig. [l] The (frequency- 
dependent) two-dimensional conductivity of graphene 
is a{uj). The ribbon can be either an actual graphene 
strip, or "virtually" created by spatially varying ex- 
ternal gates acting on a graphene sheet, as proposed 
in Ref. 12, Here we only report on the case where 
(7 = outside the ribbon, but we have checked that 
our results are virtually unmodified if this constrain is 
relaxed, provided GSPs are not supported by graphene 
for \x\ > w/2. 

Due to translational symmetry, the electric field 
of each EM eigenmode in the ribbon has the form 



Re(^),ZA 




All calculations have been performed by using the 
finite-elements commercial software "COMSOL". 



FIG. 2: (Color online) Normalized wavevector of GSP 
modes (continuous curves) and propagation lengths (dis- 
continuous curves) as a function of frequency for free- 
standing graphene ribbons of 5 /xm (a) and 20 /xm (b) 
widths. The curve for dispersion law "EGSP" is for edge 
GSPs in a semi-infinite graphene sheet, while the one 
marked as "2DGSP" is for GSP in an infinite graphene 
sheet. The colorplots for x- and z-components of the elec- 
tric field are shown in the vicinity of the ribbon cross- 
section. The numbers next to the colorplots correspond 
to the labels on the dispersion curves. 



E{f^t) = E{x, z) ex.-p{iqk^y) ex.p{—iujt)^ where k^ = 
u/c is the free space wavevector, c is the speed of 
light and q{uj) is the modal wavevector in dimension- 
less units. We will refer to the real part of q{uj) as 
the " normalized wavevector" , while its imaginary part 
provides the propagation length of GSPs, L, thorough 
L/X = 27r/Im(g'), where A = 27r/k^ is the wavelength 
in vacuum. 

For the calculations, we use the conductiv- 
ity computed within the standard random phase 
approximatioiP^H^, which depends on temperature 
(T), chemical potential (/i) and scattering time (r). 
We chose T = 300K, /i = 0.2 eV and, unless otherwise 
stated, a relaxation energy E^- = 27Th/r = O.lmeV 
(corresponding to a mobility of 1.87 x 10^ cm'^V~^s~^). 




FIG. 3: (Color online) Difference between the normalized 
wavevectors of the two coupled EGSP ribbon modes, Aq, 
together with Re(^'^^) and Im(^'^=^) (corresponding to the 
edge EM modes of a semiinfinite graphene sheet). The 
continuous, dashed and dash-dotted curves correspond to 
w = 5,10 and 20 /xm, respectively. The white circles indi- 
cate the points where Im(^'^) = A^, when the modes can 
be considered as degenerated. 

In Fig. [2] both the dispersion relations and the prop- 
agation lengths of GSPs are rendered, for two widths 
of free-standing graphene ribbons. These widths have 
been chosen so that, in the considered frequency range, 
w/X < 1 for the 5 /im- width waveguide and w/X ^ 1 
for the 20 /im-width one. Notice that the dispersion 
relation is given in the form q{i') with u = 27rcj, so 
that q = 1 sets the position of the light cone. The nor- 
malized wavevectors are shown by continuous curves, 
while the propagation lengths are represented by dis- 
continuous ones. 

We would like to notice that the microscopic struc- 
ture of graphene edges modifies significantly the elec- 
tronic spectrum, especially at energies near the Dirac 
point. The calculations reported here are not affected 
by the microscopic details of the edges. We con- 
sider graphene stripes where the separation between 
the Fermi energy and the Dirac point is much larger 
than other scales, like temperature or relaxation en- 
ergy. At these energies, the main differences between 
different types of edges is the amount of intervalley 
scattering that they induce, which vanishes for a zigzag 
edge and is maximum for armchair edges. The macro- 
scopic analysis discussed here can be written in terms 
of contributions from the two valleys. However, plas- 
mons are collective excitations of the total electronic 
charge. An electric current from a given valley incident 
at a zigzag edge is reflected in the same valley, while 
an armchair edge leads to a change of valley upon re- 
flection. The total current is, obviously, conserved in 
both cases. As the calculation depends only on to- 
tal charges and total currents, the nature of the edge 
cannot change the results. An influence of the edge 
structure can be expected, however, near the Dirac 
point, where the combination of quantum confinement 



and microscopic structure leads to changes in the elec- 
tronic spectrum. 

The dispersion relation for GSP modes in rib- 
bons can be related to that for (i) 2DGSP modes in 
an infinite graphene sheet (with wavevector q^^ = 
V"! - l/a2, where a = 27rcr/c), see curve "2DGSP" 
in Fig. ^ and (ii) surface plasmon modes appearing at 
the edges in a semi-infinite graphene sheet^^ (EGSPs) 
with wavevector q^{uj) (curve "EGSP" in Fig. 2). No- 
tice that the dispersion curve for a EGSP lies above 
that for a 2DGSP, which implies that the former is 
more tightly confined. 

We find that in a ribbon there are two modes that 
originate from the hybridization of EGSPs. Due to 
splitting, one of them (the one corresponding to even 
parity of E^ with respect to the ribbon axis), has a 
wavevector larger than q^ and is the fundamental one 
in the ribbon, while the other one has q < q^ . As seen 
from this figure, for both cases there is a single-mode 
region (shaded) for low frequencies. The wavevector of 
the fundamental mode is larger than that of 2DGSPs, 
and its propagation length is smaller, so L remains 
of the order of 10 GSP wavelengths throughout the 
considered frequency range. For sufficiently large fre- 
quencies {u = 3.8 THz for w = 5 jam and z/ = 1.8 THz 
for w = 20/im) a second GSP mode is sustained by the 
ribbon. As the frequency increases, the dispersion re- 
lations for these two GSP modes approach each other, 
merging into the dispersion relation for an EGSP of a 
semi-infinite sheet. For a given frequency, the smaller 
the width of the ribbon, the more the edge modes over- 
lap and, correspondingly, the larger the splitting be- 
tween their wavevectors q^~^ and q^~ . The highest 
frequency for which the two edge modes can be con- 
sidered as coupled can be estimated from the relation 
Aq^ < max {lm(g'^±)} c^ lm.{q^), which means that 
the "line-width" (due to absorption) of the isolated 
EM edge mode prevails the splitting (see Fig. |3|. 

Additionally to these hybridized edge modes, the 
ribbon may support several WGSPs having a field 
that extends over the ribbon width. The number of 
WGSPs trapped by the ribbon can be estimated as 
Nuj ^ 2q^^w/\^ which counts basically how many 
GSP half-wavelengths fit across the ribbon. These 
"bulk" modes present a dependence for q{uj) that lies 
below the corresponding one for a 2DGSP, which is 
their high frequency asymptote, see Fig. [2] Never- 
theless, the propagation length of these bulk modes is 
always smaller that the one for the 2DGSPs (see Fig. 

i- 

As seen from Fig. [2J for all EM modes in 
the graphene ribbon the propagation length decays 
strongly close to the modal cutoff, being much smaller 
than the one in a 2D graphene sheet, and coincides 
with the latter for high frequencies. Between these 
two regions, each curve for L/\ has a maximum at a 
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FIG. 4: (Color online) Dependencies of Re(g) and L (dis- 
continuous lines: EGSP of a semi-infinite sheet, continuous 
lines: 2DGSP) with: (a) scattering time (in units of meV), 
(b) dielectric constant of the substrate, (c) frequency. Un- 
less otherwise stated, £2 = 1 and Er — O.lmeV. 



finite frequency. To understand such a peculiar behav- 
ior let us first turn to a very interesting property of 
2DGSPs. In the studied frequency range the real part 
of the conductivity (responsible for the dissipation) de- 
creases as the frequency increases. However, due to 
increase of Re(g'), the confinement increases so quickly 
with the frequency increase that GSPs become more 
absorptive, and consequently the propagation length 
of 2DGSPs decreases. Returning to the modes in the 
ribbons, in the region of high frequencies, where their 
Re(g') are large and, actually, coincide with Re(g'^^), 
the confinement of the modes is high and their absorp- 
tion increases with the frequency, exactly as in the case 
of 2D GSPs. For lower frequencies, where Re(g') is not 
so large, the confinement is not so sensible to the value 
of q and the losses of the modes follow the character of 
the losses in a graphene sheet: they increase with the 
frequency decrease. Importantly, while q decreases, 
there are also finite wavevector components in the in- 
plane directions across the ribbon that contribute to 
the confinement. Thus, the maxima in the propaga- 
tion length corresponds to the optimum balance be- 



tween strong plasmon confinement in a less lossy sheet 
(high frequencies) and smaller confinement in a more 
absorbing sheet (lower frequencies). 

We now consider how the WGSP characteristics de- 
pend on different parameters. Their dependence on 
the relaxation energy is represented in Fig. [4] (a), show- 
ing that, while propagation length of the modes is very 
sensitive to E^ in this frequency regime, their wavevec- 
tor is practically unaffected by it. This is related to the 
fact that the imaginary part of the conductivity is al- 
most independent upon r, while the real part strongly 
depends on it, especially for lower frequencies. With 
the frequency increase, the GSPs losses become less ef- 
fected by the relaxation time of electrons being more 
sensitive to the temperature. 

With respect to the presence of a substrate, our cal- 
culations show that the dispersion relation of WGSP 
has the same structure as that shown in Fig. [2] for 
the free standing case. Figure. [4Fb) renders, as a func- 
tion of the dielectric constant of the substrate (£2), the 
spectral value of the "asymptotes", i.e. the wavevector 
and propagation length for both the EGSP supported 
by a semi-infinite graphene sheet (discontinuous lines) 
and the 2DGSP (continuous lines). As £2 increases, 
the GSPs become more localized [Re{q) increases and, 
correspondingly, the propagation length decreases]. In 



the region of high frequencies, where the momentum 
is very large, q'> l^ the GSPs are less sensitive to the 
dielectric constant of the substrate. 

To conclude, we have studied the dispersion charac- 
teristics of SP modes existing in graphene ribbons in 
the THz frequency range. We have shown that there 
are two types of SP modes: the waveguide- type, with 
the field concentrated along the whole area of the rib- 
bon (in x-direction) and the edge modes, with the field 
concentrated on the rims of the ribbons, x = ±k;/2. 
The waveguide and edge modes are separated by a 
wavevector gap from each other. The number of SP 
modes supported by the ribbon increases as either the 
frequency or the ribbon width increase. The wavevec- 
tor of the modes is not sensible to the relaxation time of 
charge carriers, but the propagation length is strongly 
affected. The SP modes can be tuned by changing 
the dielectric environment of the ribbon. The high lo- 
calization of graphene EM edge modes can be useful 
for, for instance, bending of EM signals on subwave- 
length scales and enhancing the EM coupling between 
objects. 
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